Evolution of the pygmy dipole resonance in Sn isotopes 
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Nuclear level density and y-ray strength functions of '^' '^^Sn below the neutron separation energy are ex- 
tracted with the Oslo method using the (^He.^'He'y) and Clle,ay) reactions. The level densities of '^' '^'Sn 
display step-like structures, interpreted as signatures of neutron pair breaking. An enhancement in both strength 
functions, compared to standard models for radiative strength, is observed in our measurements for Ey > 5.2 
MeV. This enhancement is compatible with pygmy resonances centered at x 8.4(1) and x 8.6(2) MeV, respec- 
tively, and with integrated strengths corresponding to x 1.8^j% of the classical Thomas-Reiche-Kuhn sum rule. 
Similar resonances were also seen in "*""'Sn. Experimental neutron-capture cross reactions are well repro- 
duced by our pygmy resonance predictions, while standard strength models are less successful. The evolution 
as a function of neutron number of the pygmy resonance in "^"'^^Sn is described as a clear increase of centroid 
energy from 8.0(1) to 8.6(2) MeV, but with no observable difference in integrated strengths. 



PACS numbers: 21.10.Ma, 24.10.Pa, 24.30.Gd, 27.60.+j 

I. INTRODUCTION 

The level density and the y-ray strength function are av- 
erage quantities describing properties of atomic nuclei. The 
nuclear level density is defined as the number of energy levels 
per unit of excitation energy, while the y-ray strength function 
may be defined as the reduced average transition probability 
as a function of y-ray energy. The strength function charac- 
terizes average electromagnetic properties of excited nuclei. 

The level density and the strength function are impor- 
tant for many aspects of fundamental and applied nuclear 
physics. They are used for the calculation of cross sections 
and neutron-capture («,y) reactions rates, which are input pa- 
rameters in, e.g., reactor physics, nuclear waste management, 
and astrophysical models describing the nucleosynthesis in 
stars. 

Tin and other heavier neutron-rich nuclei are often found 
display a smaller resonance for y-ray energies below the Gi- 
ant Electric Dipole Resonance (GEDR). The existence of even 
a small resonance close to the neutron separation energy may 
have large consequences in nuclear astrophysics on the calcu- 
lated distribution of elemental abundance. 

This article presents the measurements of the level densities 
and y-ray strength functions in '^'''-^Sn for energies below 
the neutron separation energy, as well as a systematic study of 
the evolution of the pygmy resonances in n^-i 19,121, i22(gjj 
experimental results on ' '^-"''Sn are published in Refs. HKU. 
All experiments have been peiformed at the Oslo Cyclotron 
Laboratory (OCL). 

The experimental set-up is described in Sec. II and the data 
analysis in Sec. III. The level densities of '^' '^^Sn are pre- 
sented in Sec. IV and the strength functions in Sec. V. Section 
VI discusses the pygmy resonance evolution and the impacts 
from the pygmy resonances on the (y, «) cross sections. Con- 
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elusions are drawn in Sec. VII. 



II. EXPERIMENTAL SET-UP 

The self-supporting '^^Sn target was enriched to 94% and 
had a mass thickness of 1.43 mg/cm^. For five days the 
target was exposed to a 38-MeV -^He beam with an aver- 
age current of « 0.2 nA. The reaction channels studied were 
'22Sn(3He,3He'y)'^^Sn and ^^^Sni^Hs^ajY^^Sn. 

Particle-y coincidences were recorded with 64 Si particle 
AE - E telescopes and 28 collimated Nal(Tl) y-ray detectors. 
The AE and E detector thicknesses are approximately 130 yum 
and 1550 fim, respectively. These detectors cover the angles 
of 40 - 54° with respect to the beam axis, and they have a total 
solid-angle coverage of ^ 9% out of 47r. The Nal detectors 
are distributed on a sphere and constitute the CACTUS mul- 
tidetector system ID. The detection efficiency is 15.2%, and 
the resolution of a single Nal detector is ^ 6% FWHM, at the 
y energy of 1332 keV. 



III. DATA ANALYSIS 

The measured energy of the ejectile is calculated into ex- 
citation energy of the residual nucleus. The y-ray spectra are 
unfolded with the known response functions of CACTUS and 
by the use of the Compton subtraction method |5|. The first 
generation y-ray spectra are extracted by the subtraction pro- 
cedure described in Ref. ||6|. 

The first-generation y-ray spectra are arranged in a two- 
dimensional matrix P(E, Ey), shown for '^^Sn in Fig. [l] The 
entries of the matrix give the probabilities P(E, Ey) that a y- 
ray of energy Ey is emitted from a bin of excitation energy 
E. 

The empty region for low y energy and higher excitation en- 
ergies in Fig. [T] is explained by too strong subtraction caused 
by the strongly populated states (yellow/red spots) found for 
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FIG. 1: (Color online) The first-generation matrix P of '^-Sn. The 
solid lines indicate the region for the data used in the Oslo method 
{E-y > 1.6 MeV, 3.8 < £ < 9.4 MeV). The dashed line (E = 6.6 
MeV) is the middle between the upper and lower ranges in excitation 
energy. 



low y-energy and lower excitation energies. Too few first- 
generation y's remain for low y energy and higher excita- 
tion energies, which has made the first-generation method not 
work very well (see Ref |7 1). We select and proceed with the 
region between the solid lines. It is commented that the diag- 
onal valleys and ridges are made up by strong first-generation 
transitions to the ground and first-excited states. 

The selected region of the first generation matrix P is fac- 
torized into the level density p and the radiative transmission 
coefficient 7" 10: 



P(E, Ey) OC p{E - Ey)T{Ey). 



(1) 



The factorization into two independent components is justi- 
fied for nuclear reactions leading to a compound state prior to 
a subsequent y decay |9|. The factorization is performed by 
an iterative procedure |8 | where the independent functions p 
and T are adjusted until a global minimum with the exper- 
imental P{E, Ey) is reached. 

The quality of the factorization of level density and strength 
function is illustrated for '^^Sn in Fig. |2] At example excita- 
tion energies (indicated on the panels), the entries of the P 
matrix obtained from the -fitted output functions p and T 
using Eq. ([T]) are compared to those of the experimental P ma- 
trix. The fitted output (solid line) agrees well with experimen- 
tal data. It is noted that in some of the panels, the fitted curves 
are significantly lower than the experimental data points (For 
E - A.\ MeV: the transition to the ground state; for E - 4.8 
MeV: the transition to the first-excited state). This is probably 
explained by the fit adjusting the entire matrix, and not just 



these example excitation energies. 

The Brink- Axel hypothesis llOl [m states that the GEDR 
and any other collective excitation mode built on excited states 
have the same properties as those built on the ground state. 
Equation[T]shows that the transmission coefficient is assumed 
to be independent of excitation energy E, which is a conse- 
quence of the Brink- Axel hypothesis. 

Figure [3] shows an investigation of this assumption for 
'^^Sn, which is of special concern due to some clear struc- 
tures in the strength function. We divide the selected region 
of the P matrix into two parts (separated by the dashed line 
in Fig.[T]), which are two independent data sets. Figure[3]dis- 
plays the strength functions derived from the lower and upper 
part, as well as from the total region. The strength functions, 
proportional to T/Ey^, are not normalized and are shown in 
arbitrary units. As the clear structures are found at the same 
locations for the two independent data sets, the T is indeed 
found to be approximately independent of excitation energy. 

The adjustment to Eq. ([T} determines only the functional 
forms of p and T. These two functions are invariant under the 
following transformations |[8|: 

p(E -Ey) = A exp [a [e - Ey)] p{E - Ey) , (2) 
t(Ey) = Be\p[aEy)TiEy). (3) 

The parameters A and B define the correction to the absolute 
values of the functions p and T, respectively, while the param- 
eter a defines their common correction to the log-scale slope. 
These parameters will be determined in Sees. IV and V. 



IV. LEVEL DENSITIES 

The constants A and a needed to normalize the experimen- 
tal level density p, are determined using literature values of the 
known discrete energy levels at low energy and of the level 
spacing D at the neutron separation energy 5„. We use the 
same normalization procedure as in Refs. ||T}j3|, in order to 
have a common ground for comparison. 

The chosen model is the back-shifted Fermi gas (BSFG) 
model, published by von Egidy et al. in 1988 |12|. The level 
density at the neutron separation energy p(5„) is calculated 
from known values of the ^-wave level spacing Do El: 



P(Sn) = 



Do 
+ I, exp 



(/, -H l)exp 



-a + ir 



2(r2 



2cr2 



(4) 



where /, is the target spin, and where the spin-cutoff parameter 
cr is also evaluated atS „. The spin-cutoff parameter is defined 
as o" = 0.0888 A^/^ar, where A is the mass number of the 
isotope, and T is the nuclear temperature given hyT- yfUja. 
Here, U is the nucleus intrinsic excitation energy, and a is the 
level-density parameter. The parameterization used for a is 
a - 0.21 A"**^ MeV"'. The parameterization of f/ is f/ = 
E - /ipair - C\, where the pairing energy Epair is calculated 
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FIG. 2: (Color online) Comparison between experimental (squares) and;^f^ fitted (solid lines) P matrix for '-^Sn. The energy bins have been 
compressed to 240 keV/ch in E and in Ey. While the panels show the results for the indicated example excitation energies, the fit has been 
performed globally for the entire region of the P matrix that has been selected for the analysis (see Fig.[T](. 
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FIG. 3: (Color online) Comparison of unnormalized y-ray strength 
functions (arbitrary units) for '-^Sn. The shown strength functions 
are derived from the independent data sets of the upper and lower 
part of the selected region of the P matrix, as well as from the total 
selected region for comparison. 



from the proton and neutron pair-gap parameters: Ep^ii- - Ap 
An, and where the back-shift parameter Ci is defined as Ci 



The experimental value of Dq for Sn is found in Ref. |fT3]| 
and is used to calculate p{S „) using the input parameters 
listed in Tab. [l] The pair-gap parameters are evaluated from 
even-odd mass differences lfT4l according to the method of 
Ref. [T5l. 

No experimental value exists for Dq of '^^Sn, and we esti- 
mate p{S„) for this isotope from systematics. Figure |4] shows 
p{S „) calculated from the experimental values of Do accord- 
ing to Eq. (|4]) for all available Sn isotopes as a function of S„. 
The values of Do have been taken from Ref. |13|. We have 
also calculated p{S„) according to the prediction of the BSFG 
model I.16J : 



P(£')bsfg 



exp (2 yfaU) 
12V2fli/4 t/5/4o- 



(5) 



with the above-listed parameterizations. The theoretical value 
of p(S„), multiplied with a common factor of 0.4, are shown in 
Fig. [4] together with the experimental values. From the trends 
appearing in the figure, we estimate p{S„) for '^^Sn to 2.0(10)- 
10^ MeV"' (50% uncertainty assumed, see Tab.|l]i. 

While we would like to normalize to p(S„), our experimen- 
tal data only cover the excitation energy region from to 5 „ -2 
MeV, due to methodical limitations. We therefore make an 
interpolation from our measurements to S„ using the BSFG 
prediction in Eq. (|5]l. The prediction is multiplied by a scaling 
parameter 77 (see Tab. in order to agree with the normaliza- 
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TABLE I: Input parameters and the resulting values for the calculation of the normalization value p(5„). 



Nucleus 


5„ Oo a C| A„ Ap cr(5„) p(5„) 
(MeV) (eV) (MeV"') (MeV) (MeV) (MeV) (10* MeV"') 


1 


"'Sn 
'22Sn 


6.17 1250(200) 13.62 -1.42 0.82 4.57 3.42(86) 
8.81 62(31)° 13.72 -1.42 1.32 1.12 4.75 20(10)° 


0.25 
0.42 



Estimated from systematics. 
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FIG. 4: (Color online) Estimation (cross) of the experimental value 
of p{S „) for '-^Sn from systematics. Experimental values (squares) 
and global BSFG predictions (triangles) for p(5„) are shown for var- 
ious Sn isotopes as function of 5„. (See text.) 



tion value p(S„y. 
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(6) 



Figure [5] shows the normalized level densities of '^' '^^Sn. 
The arrows indicate the two regions that have been used for 
normalization to the discrete level density and to the normal- 
ization value piS„). As expected, the level density function 
of '^'Sn is very similar to those of the other even-odd nuclei 
ii7,ii9(5jj^ while the level density function of '^^Sn is very sim- 



show step-like 
1. In I2i.i22sn, 



ilar to those of the even-even nuclei ''^ "^Sn 
The level densities of '^^ '^^Sn in Fig. p 
structures, a feature also seen in "^""''Sn im 
two pronounced bumps are seen in the region of ^ 0.8 - 1.4 
MeV and a; 1.8 - 2.4 MeV. The corresponding steps are lo- 
cated at a; 1.0 and x 2.0 MeV, respectively. The second step is 
very abrupt, especially in the even-even nucleus, and the step 
is followed by a significantly higher level density. The sec- 
ond step is therefore a candidate for the neutron pair-breaking 
process in '^'-'^^Sn. Such neutron pair-breaking bumps are 
especially distinctive in Sn isotopes since the proton number 



10" 



10=" 



Si 3 
>. 

1 10^ 

J 10 



10=" 



iio^ 



10" 



Q. 



£10^ 
V 
T5 

1 10 



10-^ 



1 r 



121 



Sn 



I I ... I .... I .... I .. .1. 1 1 ... I .... I 



122 



Sn 



.---'1 



ni^ i 



2 3 4 5 6 7 

Excitation energy E (lUleV) 



FIG. 5: Normalized level densities of '-'-'^"Sn (filled squares) as a 
function of excitation energy with energy bins of 120 keV/ch. The 
solid lines represent the discrete level densities obtained from spec- 
troscopy. The dashed line in both panels is the BSFG prediction, 
which is used for interpolation, scaled with rj to coincide with p(S„) 
(open square). The value of p(5„) has been calculated from neutron 
resonance data. The arrows indicate the two regions used for nor- 
malization. 



is magic (Z = 50), making proton pair breaking occur only at 
relatively higher excitation energies. A detailed discussion of 
the pair-breaking process has been given in Refs. |[T][3l. 
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V. GAMMA-RAY STRENGTH FUNCTIONS 



The y-ray transmission coefficient 7", which is deduced 
from the experimental data, is related to the y-ray strength 
function / by 



T{Ey)^2nyEl^^\fxL{Ey), 



XL 



(7) 



where X denotes the electromagnetic character and L the mul- 
tipolarity of the y ray. The transmission coefficient T is nor- 
malized in log-scale slope (a) and in absolute value {E) (see 
Eq. 0). 

For i-wave neutron resonances and assuming a major con- 
tribution from dipole radiation and parity symmetry for all 
excitation energies, the expression for the average radiative 
width (Ty{E, I, n)^ will at 5„ reduce to lil7J : 

{Yy(S„J,±\l2,n,)) 

= A re f^Mi ^ f " dEyT(Ey)p(S„-Ey) 
47rp(5„,/f + l/2,7r,) Jo 
1 

xY,8(Sn-Ey,I,±l/2 + J). (8) 



Here, /, and n, are the spin and parity of the target nucleus in 
the neutron capture in,y) reaction. We determine B by using 
the BSFG model for the spin distribution g given in Ref. Iil2l 
and the experimental value of (ry{S„)y 

For '^'Sn, the radiative width at the neutron separation en- 
ergy is available in literature. For '^^Sn, we estimated it 
from systematics. Figure 6 shows the (ry(S„)^ plotted against 
S„ for Sn isotopes wherethis quantity is known (taken from 
Ref. fT3l). From the appearing trend of the even-even nuclei, 
we estimate (^ry{S„)^ to 85(42) meV for '^^Sn. The applied in- 
put parameters needed for determining the normalization con- 
stant B for '^'"'^^Sn are shown in Tab.Hll The values for '^'Sn 
have been taken from Ref. ifTSll . 



TABLE II: Input parameters for normalization of the y-ray transmis- 
sion coefficient T of '^' '^^Sn. 



Nucleus 


/, 


Do 


{TyiS,,)) 




(%) 


(eV) 


(meV) 


"'Sn 





1250(200) 


40(5) 


'22Sn 


3/2 


62(31)" 


85(42)" 



160 

I 140 

c 

(A 

7a 120 

£ 100 

•a 



> 

(0 
(0 

(0 



a> 40 



80 



60 
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X Estimated 
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FIG. 6: Estimation of the average radiative width at 5„, {T-yiS „)j, 
for '^^Sn. The respective values of other Sn isotopes are plotted 
(squares) as a function of 5„. (See text.) 



" Estimated from systematics. 



onset of a small resonance, commonly related to as the pygmy 
dipole resonance. A comparison of our '^'''-^Sn measure- 
ments compared with photonuclear cross section data from 
Refs. fTS'-STi is shown in the two upper panels in Fig. [8] Sim- 
ilar strength increases were also seen in '""""^Sn |3 1, and this 
figure will be further discussed for those isotopes when dis- 
cussing the evolution of the pygmy resonance in the next sec- 
tion. 



In order to investigate the experimental strength functions 
of '^' '^-Sn, we have applied commonly used models for 
the GEDR resonance and for the magnetic spin-flip reso- 
nance, also known as the Giant Magnetic Dipole Resonance 
(GMDR). 



The normalized y-ray strength functions of '^' '^^Sn are 
shown in Fig. |7] The error bars show the statistical uncer- 
tainties. While the strength function of '^'Sn is smooth, just 
like those of "''Sn EH, the strength function of '^^Sn 
displays clear structures in the entire Ey region. As discussed 
in Sec. Ill, these structures also appear using different, inde- 
pendent data sets. 

A change of the log-scale slope in the strength functions, 
leading to a sudden increase of strength, is seen in '^' '^^Sn for 
Ey > 5.2 MeV. The change of log-scale slope represents the 



The Generalized Lorentzian (GLO) model |23 | is used for 
the GEDR resonance. The GLO model is known to agree 
rather well both for low y-ray energies and for the GEDR cen- 
troid at about 15 MeV. The strength function approaching a 
non-zero value for low Ey is not a property specific for the Sn 
isotopes, but has been seen for all nuclei studied at the OCL 
so far 



In the GLO model, the £1 strength function is given by 
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FIG. 7: Normalized y-ray strength functions of '-'"'^-Sn as functions 
of y-ray energy. The energy bins are 120 keV/ch. 



^KMpiEy, Tf) 



+ 0.7 



^KMF(Ey = 0, Tf) 



(9) 



in units of MeV"^, where the Lorentzian parameters are the 
GEDR's centroid energy Ee\, width F^i and cross section 
(Tei- These experimental parameters are not available for 
I2i,i22(gjj instead apply the one measured for '^"Sn to 
'^'Sn, and the one measured for '^"^Sn to '^^Sn, from Fultz 
ifTSl (seeTab.|IlI|. 

The GLO model is temperature dependent from the incor- 
poration of a temperature dependent width F^mf- This width 
is the term responsible for the non-vanishing £1 strength at 



TABLE III: Applied parameters for the parameterization of the 
GEDR and the GMDR contributions for '"' '^-Sn. 



Nucleus Ef 



o-Ei 



O" Ml 



(MeV) (MeV) (mb) (MeV) (MeV) (mb) (MeV) 



'^'Sn 
'22Sn 



15.53 
15.59 



4.81 
4.77 



253.0 
256.0 



8.29 
8.27 



4.00 
4.00 



1.11 0.25(5) 
1.09 0.25(5) 



low excitation energy. It has been adopted from the Kadmen- 
skii, Markushev and Furman (KMF) model I.24J and is given 
by: 

TKMFiEy, Tf) = ^ [Ey^ + An'Tf'-) , (10) 

in units of MeV, and where Tf is the temperature. 

Usually, Tf is interpreted as the nuclear temperature 
of the final state, with the commonly appUed expression 
Tf - yfUja. In this work and in Refs. lU-O, we assume 
a constant temperature, i.e., the y-ray strength function is in- 
dependent of excitation energy. This approach is adopted for 
consistency with the Brink- Axel hypothesis (see Sec.[III|. 

Moreover, we treat Tf as a free parameter This is necessary 
to adjust the theoretical strength prediction to our low-energy 
measurements. The applied values of the Tf parameters are 
listed inTab.lm] 

The Ml spin-flip resonance is modeled with the functional 
form of a Standard Lorentzian (SLO) model I 25J : 



1 



Ml 



Ey - Em\ ) + Ey Tm\" 



(11) 



where the parameter Em\ is the centroid energy, F^i the width 
and ctmx the cross section of the GMDR. These Lorentzian 
parameters are for '^' '^^Sn predicted from the theoretical ex- 
pressions in Ref. [25 1 and shown in Tab. Ill The predictions 



for the GEDR using the GLO model and for the GMDR for 
116-119,121, i22jgjj nuclei are shown as dashed lines in Fig.js] 

The Standard Lorentzian (SLO) model was also tested and 
is included in Fig.[8](the Ml spin-flip resonance contribution 
is added to it). The SLO succeeds in reproducing the (y, x) 
data, but clearly fails for the low-energy strength measure- 
ments, both when it comes to absolute value and shape. The 
same has been the case also for many other nuclei measured 
at the OCL and elsewhere. Therefore, we consider the SLO 
not to be adequate below the neutron threshold. 

At present, it is unclear how these resonances should be 
modeled properly, although many theoretical predictions ex- 
ist. We have found i21|3J that the Sn pygmy resonance is sat- 
isfactorily reproduced by a Gaussian distribution [2J: 



fpyg(.Ey) — C. 



1 



pyg 



^/27T^ 



exp 



• (^Ey ^pyg) 

2a- 2~ 



2\ 



(12) 



superimposed on the GLO prediction. Here, Cpyg is the res- 
onance's absolute value normalization constant, Epyg the cen- 
troid energy and o-^yg the width. These parameters are treated 
as free. 
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FIG. 8: (Color online) Comparison of the prediction of the total strength functions with the OCL experimental measurements for 
116-119,121. i22gjj jj^g strength function of "^Sn is shown both as multiplied with 1.8 (filled squares) and as the original normalization (open 
squares). The arrows indicate S„. The total strength predictions (solid lines) are modeled as Gaussian pygmy resonance additions to the GLO 
(El + Ml) baselines. The SLO (El + Ml) baselines are also shown. 

Upper, left panel: ' ' ''Sn(7, n) from Utsunomiya et al. 1 19 1, ' "■'^'*Sn(7, x) from Fultz et al. 1181 , and '^''Sn(7, x) from Varlamov et al. 1211 . 
Upper, right panel: "*Sn(7,«) |T9l, '^''•'2^Sn(y, x) 1 18|, and ^^^S\i{y,n) from Varlamov et al. (jp). 
Middle, left panel: "'Sn(7,n) (T9|, "''Sn(7,A:) QlJ, and "''Sn(7,n) 120|. 

Middle, right panel: '"'Sn(7, n) |T9|, '"*Sn(7, x) (HI, '"*Sn(7, x) from Varlamov et al. (21], and ' '**Sn(7, x) from Lepretre et al. l22l . 
Lower, left panel: '"Sn(7,n) (191, '"Sn(7,x) (HI, '"Sn(7,x) (21], and '"Sn(7,x) l22l . 
Lower, right panel: "'^Sn(7,w) (19], '"'Sn(7,x) (H], '"'Sn(7,x) O, and '"^Sn(7,x) (22]. 



By adding the discussed theoretical strength contributions, 
the prediction of the total y-ray strength function is given by: 



/t 



total 



(13) 



fsi + /mi + /pyg • 

We determined the Gaussian pygmy parameters of '^'"'^^Sn 
from fitting to our measurements. The centroid energies of the 
pygmy resonances are 8.4(1) and 8.6(2) MeV, respectively. 
We found that the same width o-pyg and strength Cpyg as in 
ii6,ii7(gjj |3| gave a very good agreement also in ''^^Sn, 
so the width and strength are kept unchanged. The pygmy 



parameters are listed in Tab. |IV] The estimated error bars 
given in the table take into account systematic uncertainties in 
the normalization values and in the choice of baseline of the 
pygmy resonance, including the fact that the GLO does not 
perfectly follow the (y, n) measurements for higher Ey values. 

The predictions for '^' '^-Sn are shown as solid lines in the 
upper panels of Fig.|8] We see that the predictions of the total 
strength describe the measurements rather well, in the sense 
that the Gaussian pygmy resonances fill a very large fraction 
of the missing strength. Still, the Gaussian distribution does 
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TABLE IV: Empirical values of Gaussian pygmy parameters, and 
the corresponding integrated strengths and TRK values of the pygmy 
resonances, in '^' '^^Sn. 



Nucleus E„. 



C 



Integrated strength TRK value 



(MeV) (MeV) (lO^'MeV"^) (MeV-mb) 



'^'Sn 
'22Sn 



8.4(1) 
8.6(2) 



1.4(1) 
1.4(1) 



3:2!^ 

3.2;3 



31+y 
32!^ 



(%) 



+T 

1.8+1 



not completely follow neither the left flank nor the right flank 
of flie pygmy resonances in '^'"'^^Sn. In the case of "^'"^Sn 
Il2l|3l, having a larger Tj of Tf = 0.46(1) MeV, the left flank 
was followed very well. However, in all Sn isotopes, there is a 
gap on the right flank between measured data and the GLO. A 
better pygmy resonance representation than the Gaussian or a 
better model for the baseline may be found in the future. 

Strength from the resonances in '^' '-^Sn have been added 
in the energy region of » 5-8 MeV according to our measure- 
ments, and in the region of ^ 5-11 MeV when comparing 
to photonuclear data as well. The total integrated strengths 
of the pygmy resonances, based on the total predictions, are 
estimated to 32'^l MeV-mb. This constitutes 1.8+5% °f 
classical Thomas-Reiche-Kuhn (TRK) sum rule, assuming all 
pygmy resonance strength being £1, see Tab.|IV| 

Even fliough uncertainties in the choice of baseline have 
been considered in the uncertainty estimates, another predic- 
tion of the GEDR than the GLO or another function for the 
pygmy resonance than the Gaussian may be found in the fu- 
ture. This will consequently influence the estimates on the 
pygmy resonance. Lack of data, i.e., the gaps between our 
measurements and the (y, n) measurements in resonance re- 
gion, and also the lack of (y, n) measurements for '^'"'^^Sn, 
adds to the uncertainties in the estimates of the pygmy reso- 
nances. 

Measurements from other reactions and using other meth- 
ods have also been used to estimate the TRK value of the Sn 
pygmy, and these estimates deviate from each other. Data 
from "^'"^Sn(y,n) experiments |19| indicate an integrated 
strength which constitutes » 1% of the TRK sum rule, 
which agrees within the uncertainty with our value. From 
'"'•'^'*Sn(y,7'') experiments [26], the TRK value is calculated 
to 0.4 - 0.6%. This may seem to deviate from our result. 
However, taking into account unresolved strength in the quasi- 
continuum of typically a factor of 2 - 3, the (7,7') results are 
compatible within the uncertainty with the other data. 



VI. EVOLUTION OF THE PYGMY RESONANCE 

Studying the neutron dependency is important and may 
help in determining the origin of the Sn pygmy resonance. 
Figure |9] shows the present and previously analyzed normal- 
ized strength functions for the Sn isotopes. The measurements 
of "^Sn have been multiplied by 1.8 in order to agree with 
those of ii^'Sn (see Ref |3|). 

First, it may seem like a trend that the tail of the strength 
function decreases in strength as the neutron number in- 




2 3 4 5 6 7 

Y-ray energy E (MeV) 



FIG. 9: (Color online) Normalized y-ray strength functions as func- 
tions of y-ray energy for the Sn isotopes measured at the OCL. The 
even-odd isotopes, "' "''•'^'Sn, are shown in the upper panel, while 
the even-even isotopes, "* '"* '^-Sn, are shown in the lower. The mea- 
sured strength of "*Sn has been multiplied by 1.8 (see Ref. [3 |). The 
energy bins are 120 keV/ch for 116.117.121.1223^^ 240 keV/ch for '"Sn 
and 360 keV/ch for "*Sn. 



creases. Second, it is noticeable from Fig. [9] that the change 
of log-scale slope, which represents the onset of the pygmy 
resonance, occurs at an higher Ey value in '^' '^^Sn than for 
116,1 ivi^jj -pjjg changes of slope are clearest for the even-even 
nuclei. They are found at » 4.5 MeV in "^Sn and at x 5.2 
MeV in i^^Sn. 

The values of Tf for '^' '^^Sn are lower than for "^""^Sn. 
There is no physical reason for different nuclear temperatures. 
Lowering the parameter T/ is instead necessary in order to get 
the lowest-energy part of the GLO comparable in magnitude 
with the measurements. 

The centroid energy Epy„ of the pygmy resonances in 
I2i,i22(gjj jj^g larger values than those of earlier studies in 
"^■"^Sn EH. For ^^^'^^^Sn, the pygmy centroids ire 8.4(1) 
and 8.6(2) MeV (see Tab.|IV]l, respectively, while 8.0(1) MeV 
is found for "^ "^Sn |2 3|. During the pygmy resonance fit- 
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TABLE V: Empirical values of "^""'Sn Gaussian pygmy parame- 
ters, and the corresponding pygmies' integrated strengths and TRK 
values. For "^'Sn, the values have been found from fitting to the 
measured strength function multiplied by 1.8. 



Nucleus 


F rr 

-^pyg ^ pyg 


c 


Integrated strength TRK value 




(MeV) (MeV) (10 


-^MeV-2) 


(MeV-mb) 


(%) 


"^Sn 


8.0(1) 1.4(1) 


3.2!^ 


30^8 


1.7!;; 


'"Sn 


8.0(1) 1.4(1) 


3.2!^ 


30+° 


1.7!° 


'"*Sn 


8.2(1) 1.4(1) 


3.2!° 


30+° 


1.8+° 


""Sn 


8.2(1) 1.4(1) 


3.2!° 


30!° 


1.7!: 



TABLE VL Applied parameters for the parameterization of the 
GEDR and the GMDR contributions for ' "Sn. 



Nucleus 


Eei 


Ffi 


0"£1 




Tmi 


C Ml Tf 




(MeV) (MeV) 


(mb) 


(MeV) (MeV) 


(mb) (MeV) 


I'^Sn 


15.68 


4.19 


266.0 


8.41 


4.00 


0.773 0.46(1) 


"'Sn 


15.66 


5.02 


254.0 


8.38 


4.00 


1.04 0.46(1) 


'"*Sn 


15.59 


4.77 


256.0 


8.36 


4.00 


0.956 0.40(1) 


"'Sn 


15.53 


4.81 


253.0 


8.34 


4.00 


0.963 0.40(1) 



ting, it was clear that the centroid energies had to be signifi- 
cantly increased for the heavier isotopes. The significant in- 
creases are also apparent from studying the energies for which 
there is a change of log-scale slope in the strength functions. 
Moreover, keeping the centroid energy constant has as a con- 
sequence that the same pygmy width cTpyg and pygmy strength 
Cpyg as in '"'•"'Sn |3J also give the best fit in ^'^^'^'^-Sn. 

In the earlier study of '^^''^'Sn ||3], the data have large eiTor 
bars. This means that the pattern of an increasing centroid en- 
ergy was not so clear, and the choice then was to keep the cen- 
troid energy constant while compensating with an increase of 
the resonance width. We have updated the resonance predic- 
tion of "^'"^Sn by following the same pattern. The estimated 
centroid energy of the pygmies in "^ "''Sn is then 8.2(1) MeV, 
while the width and strength are kept constant. Updated pa- 
rameter values are listed in Tab. [V] and displayed in Fig. [8] 
The parameters for the GEDR and GMDR contributions are 
listed inTab.|Vll 

We would like to investigate for several isotopes the ef- 
fect of our predicted pygmy resonances on the («,y) cross 
sections and compare these with available experimental mea- 
surements. This has been done for ii''-ii9,i2ijgjj using the re- 
action code TALYS |27|. For the level density, we have ap- 
plied the spin- and parity-dependent calculations of Goriely, 
Hilaire and Koning |28|, which are in good agreement with 
our level-density data, as demonstrated for "^Sn in Fig 
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Also, we have used the neutron optical potential of Koning 
and Delaroche f29\. 

The results of the comparisons are shown in Fig. 1 1 Our 
modeled strength function with a Gaussian pygmy resonance 
(denoted GL02) leads to a calculated cross section that gen- 
erally agrees very well with the measurements. Assuming the 
GLO model with constant temperature but without the pygmy 
resonance (GLOl), clearly gives a lower cross section in all 
cases, as expected. This may be taken as a support of the 
finding of an enhanced strength function in the Ey region of 
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FIG. 10: (Color online) Level-density measurements on "'Sn com- 
pared with microscopic calculations from Ref. 1281 . 



a; 5 - 1 1 MeV. The SLO model gives an overall too high cross 
section, which is not surprising considering the large over- 
shoot in y-ray strength compared to our measurements and 
also to ^Fy^ data. We note that our calculated cross section 

for '"'Sn(n,y)"'Sn using the GL02 model is in very good 
agreement with the one in the work of Utsunomiya et al. |19| . 

For the "'Sn(«, y)"**Sn case, we have applied the model 
parameters that coiTespond to our scaled data (with a factor of 
1.8). The resulting excellent agreement with the experimen- 
tal in, y) data further supports our choice of renormalization 
in Ref. lO. We have in addition tested the strength predic- 
tion using parameters that fit with the original normalization, 
which results in calculated cross section that are clearly un- 
derpredictive compared to the experimental data (not shown 
here). 

The reproduction of the '^''Sn(n, 7)'^'Sn cross section is not 
as good as for the other nuclei, as this calculation seems to 
be somewhat more underpredictive. However, the calculation 
is certainly an improvement compared to that of the GLOl, 
which is a standard strength model without the pygmy reso- 
nance. The underprediction might be explained by too low ex- 
perimental value of ^Fy^ in the normalization procedure of the 
measurements. If the value had been higher, the pygmy reso- 
nance would had produced more strength, leading to a general 
increase of the calculated cross section. 

We would also like to study the evolution of the resonances' 
centroid energy E^^„ with neutron number A^. 



Figure 



12 shows 



/ipyg as a function of for the isotopes studied at theOCL. A 
X fit has been performed on these data, resulting in the linear 
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FIG. 11: (Color online) Data on neutron-capture cross sections for the target nuclei ii6-ii8,i20gj^ compared with calculations, for neutron 
energies > 20 keV. GLOl (dotted line) is the GLO model assuming constant temperatures (given in Tabs. |IV| and |VI^ , and GL02 (solid line) 
is the GLOl model plus the prediction of the pygmy resonance. Measurements from Koehler et al. L30J . Macklin et at. L31J , Nishiyama et 
al. 1321 . Wisshak et al. 1331 . and Timokhov et al. 1341 . 



relation E^yg = 2.0(16) + 0.090(23) ■ in units of MeV. 

The estimates on Epyg from others' experimental data on Sn 
are in agreement within the uncertainties with the observed 
pattern: * 8.5 MeV for '"'•"^Sn 1 19|, ~ 7.8 MeV for "^•"''Sn 
Il35l . and 10.1(7) and 9.8(7) MeV for '^^ '^ZSn respec- 
tively. It is commented that an increase of the resonances' 
centroid energies with increasing neutron number was also 
found in experimental data on Ca BTI . 

The observation of an increase of the centroid energy 
with increasing neutron number is not in agreement with 
common theoretical predictions. On the contrary, studies 
on Sn isotopes instead predict a decrease of centroid en- 
ergy with increasing neutron number. These studies in- 
clude the Hartree-Fock-Bogoljubov (HFB) and multiphonon 
quasiparticle-phonon (QPM) models by Tsoneva and Lenske 
p8|, the RHBh-RQRPA (DD-ME2) model by Paar |39|, and 
the Continuum QRPA model by Daoutidis BDI . Also a the- 



oretical study on Ca isotopes, using the Extended Theory of 
Finite Fermi Systems (ETFFS) by Tertychny ll37ll . results in a 
decrease of centroid energy with neutron number (in contrary 
to experimental results on Ca, see Ref. 1371 and references 
therein). However, it is commented that Daoutidis |40| pre- 
dicts a relatively stable centroid energy in the atomic mass re- 
gion A - 120 - 126 compared to other mass regions. Hence, 
the increase of centroid energy in the isotopes that we have 
compared, may be less than would had been observed in an- 
other mass region. 

Recent measurements using the {a,a'y) coincidence 
method on '^"^Sn compared to photon-scattering experiments 
show a splitting into its isoscalar and isovector components 
||4T1| . Hence, both components seem to be present, in agree- 
ment with theoretical predictions. 

The nature, origin, and integrated strength of the Sn pygmy 
resonance are issues that are heavily debated. The El neutron- 
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FIG. 12: (Color online) Estimated centroid energies Epy, (squares) as 
a function of neutron number A', deduced from our measurements on 
116-1 19,121, 1225JJ jjjg sQjijj ^[jg is a linear fit to the measurements. 



FIG. 13: (Color online) TRK values for Sn estimated from our mea- 
surements (squares) compared to the theoretical prediction from Van 
Isacker et al. |42 1 (multiplied by a factor 14) (solid line) as a function 
of neutron number A'. 



skin oscillation mode, discussed in Refs. Il39l l42l |431 . is as- 
sumed as the underlying physical phenomenon in most of the 
theoretical predictions, both in macroscopic (e.g.. Van Isacker 
et al. |42 |) and microscopic approaches (e.g., Daoutidis BOl . 
Tsoneva and Lenske |381, Paar llU, Sarchi et al. El). Most 
theoretical calculations predict a systematic increase of the 
resonances' strength as the number of neutrons increase, due 
to the increase of the number of neutrons in the skin. An- 
other prediction is the increase by neutron number up to '^°Sn 
followed by a decrease (Paar |[39l ). Several of the predicted 
increases of integrated strength concerning the isotopes that 
we have performed measurements on, are significant (e.g., 
Tsoneva and Lenske Il38l . Van Isacker et al. Il42l . Litvinova 
et al. pTi). It is commented that the study by Daoutidis |l40l 
predicts that also the integrated strength is relatively stable in 
the mass region A = 120 - 126. 

However for our measurements on the pygmy resonances 
jjj 116-119,121, i22(gjj^ cannot see any dependency on neutron 
number in the integrated strengths. The same resonance pre- 
diction has on the contrary been applied for all isotopes. (The 
total integrated strengths and the TRK values of '^' '^^Sn be 
ing slightly larger than those of '"'•"^Sn, see Tabs 



IV 



and 



IV 



is explained by differences in the GLO models of those iso- 
topes.) Figure [T3] shows our TRK values together with those 
of Van Isacker et al. (42] (multiplied by a factor of 14 in abso- 
lute value). The experimental result does not follow the pre- 
dicted increase. Still, the uncertainties in our estimated reso- 
nance strengths are large. More experimental information is 
therefore needed in order to answer the question if the inte- 
grated strength in Sn increases with neutron number 



The experimental TRK values based on our and others' 
measurements are relatively large, compared to general ex- 
citations. This might indicate that the pygmy resonance is 
due to a collective phenomenon. However, its origin is un- 
known, and single-particle excitations are not excluded. Var- 
ious theoretical predictions disagree on whether the neutron 
skin-oscillation is collective or not \39]. 

A clarification of the electromagnetic character of the 
pygmy resonance in Sn would be of utmost importance. 
Present theoretical predictions assumes an £1 strength, mod- 
eling the resonance as a neutron-skin oscillation. Experi- 
mental studies have indicated an £1 character. Amongst 
these are the nuclear resonance fluorescence (NRF) experi- 
ments performed on "^•'24sn ||26l and on "2,124sjj gjj^ j^^^j 
the Coulomb dissociation experiments performed 
Il36l l4&il . In addition comes the polarized photon beams ex- 
periments on '-^^Ba BtI . However, the existence of an Ml 
component of the resonance strength cannot be ruled out, as 
was discussed in Ref. |i3J. 



VII. CONCLUSIONS 

The level density and y-ray strength functions of ^^' '^^Sn 
have been measured using the (^He,^He''x) and (^He,ay) re- 
actions and the Oslo method. The level densities of ^^' '^^Sn 
display step-like structures for excitation energies below « 4 
MeV One of the bumps is interpreted as a signature of neutron 
pair breaking, in accordance with the findings in "''Sn. 
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A significant enhancement in the y-ray strength is observed 
in the '^'"'^^Sn measurements for Ey > 5.2 MeV. The inte- 
grated strength of the resonances correspond to x 1.8^j% of 
the TRK sum rule. These enhancements are compatible with 
pygmy resonances centered at « 8.4(1) and » 8.6(2) MeV, 
respectively. 

Neutron-capture cross-section calculations using our 
pygmy resonance predictions give significantly better repro- 
duction of experimental («, y) cross sections than standard 
strength models without any pygmy resonance. 

The pygmy resonances are compared to those observed in 
ii6-ii9§jj -j^g evolution with increasing neutron number of 
the pygmy resonances observed in Ii6-ii9,i2i,i22|gj^ ^ clear 
increase of centroid energy from 8.0(1) MeV in ^'^Sn to 8.6(2) 
MeV in i22sn, while no difference in integrated strength is 
observed. This finding is not in agreement with most theo- 
retical predictions. However, the experimental uncertainties 



are large, and more experimental information is needed in or- 
der to determine the nature of the pygmy resonances in the Sn 
isotopes. 



Acknowledgments 

The authors would like to thank H. Utsunomiya for giv- 
ing us access to yet unpublished experimental results from 
photoneutron cross-section reactions '^°'^^^Sn(7,n). The De- 
partment of Physics at the University of Jyvaskyla (JYFL) is 
thanked for kindly lending us the '^^Sn target. Last but not 
least, thanks to E. A. Olsen, J. Wikne, and A. Semchenkov 
for excellent experimental conditions. The funding of this re- 
search from The Research Council of Norway (Norges forskn- 
ingsrad) is gratefully acknowledged. 



[1] U. Agvaanluvsan, A. C. Larsen, M. Guttormsen, R. Chankova, 
G. E. Mitchell, A. Schiller, S. Siem, and A. Voinov, Phys. Rev. 
C 79, 014320 (2009). 

[2] U. Agvaanluvsan, A. C. Larsen, R. Chankova, M. Guttormsen, 

G. E. Mitchell, A. Schiller, S. Siem, and A. Voinov, Phys. Rev. 
Lett. 102, 162504 (2009). 

[3] H. K. Toft, A. C. Larsen, U. Agvaanluvsan, A. Burger, M. Gut- 
tormsen, G. E. Mitchell, H. T. Nyhus, A. Schiller, S. Siem, N. 
U. H. Syed, A. Voinov, Phys. Rev. C 81, 064311 (2010). 

[4] M. Guttormsen, A. Atac, G. L0vh0iden, S. Messelt, T. Rams0y, 
J. Rekstad, T. F. Thorsteinsen, T. S. Tveter, and Z. Zelazny, 
Phys. Scr. T 32, 54 (1990). 

[5] M. Guttormsen, T. S. Tveter, L. Bergholt, F. Ingebretsen, and 
J. Rekstad, Nucl. Instrum. Methods Phys. Res. A 374, 371 
(1996). 

[6] M. Guttormsen, T. Rams0y, and J. Rekstad, Nucl. Instrum. 

Methods Phys. Res. A 255, 518 (1987). 
[7] A. C. Larsen, M. Guttormsen, E. Betak, A. Biirger, A. Gor- 

gen, M. Krticka, H. T. Nyhus, J. Rekstad, A. Schiller, S. Siem, 

H. K. Toft, G. M. Tveten, A. Voinov, and K. Wikan, A review of 
the Oslo method, submitted to Phys. Rev. C. (2010). 

[8] A. Schiller, L. Bergholt, M. Guttormsen, E. Melby, J. Rekstad, 
and S. Siem, Nucl. Instrum. Methods Phys. Res. A 447, 498 
(2000). 

[9] A. Bohr and B. Mottelson, Nuclear Structure (Benjamin, New 
York, 1969), Vol. I. 
[10] D. M. Brink, Ph.D. thesis, Oxford University, 1955. 
[11] P Axel, Phys. Rev 126, 671 (1962). 

[12] T. von Egidy, H. H. Schmidt, and A. N. Behkami, Nucl. Phys. 
A481, 189 (1988). 

[13] R. Capote et al, RIPL-3 - Reference Input Parameter Library 
for Calculation of Nuclear Reactions and Nuclear Data Evalu- 
ations, Nucl. Data Sheets 110, 3107 (2009). Available online at 
[http://www-nds.iaea.org/RIPL-3/]. 

[14] G. Audi and A. H. Wapstra, Nucl. Phys. A595, 409 (1995). 

[15] J. Dobaczewski, P. Magierski, W. Nazarewicz, W. Satula, and 
Z. Szymaliski, Phys. Rev. C 63, 024308 (2001). 

[16] A. Gilbert and A. G. W. Cameron, Can. J. Phys. 43, 1446 
(1965). 

[17] A. Voinov, M. Guttormsen, E. Melby, J. Rekstad, A. Schiller, 
and S. Siem, Phys. Rev. C 63, 044313 (2001). 



[18] S. C. Fultz, B. L. Berman, J. T. Caldwell, R. L. Bramblett, and 

M. A. Kelly Phys. Rev 186, 1255 (1969). 
[19] H. Utsunomiya, S. Goriely, M. Kamata, T. Kondo, O. Itoh, H. 

Akimune, T. Yamagata, H. Toyokawa, Y.-W. Lui, S. Hilaire, 

and A. J. Koning, Phys. Rev C 80, 055806 (2009). 
[20] V. V. Varlamov, B. S. Ishkhanov, V. N. OrUn, V. A. Tch- 

etvertkova, Moscow State Univ. Inst, of Nucl. Phys. Reports 

No. 2009, p. 3/847 (2009). 
[21] V. V. Varlamov, N. N. Peskov, D. S. Rudenko, and M. E. 

Stepanov, Vop. At. Nauki i Tekhn., Ser. Yademye Konstanty 

1-2 (2003). 

[22] A. Lepretre, H. BeU, R. Bergere, P. Carlos, A. De Miniac, A. 
Veyssiere, and K. Kembach, Nucl. Phys. A 219, 39 (1974). 

[23] J. Kopecky and R. E. Chricn. Nucl. Phys. A 468, 285 (1987). 

[24] S. G. Kadmenskii, V. P. Markushev, and V. I. Furman, Yad. Fiz. 
37, 277 (1983) [Sov J. Nucl Phys. 37, 165 (1983)]. 

[25] T. Belgya et ah. Handbook for calculations of nuclear reac- 
tion data, RIPL-2, (IAEA, Vienna, 2006). Available online at 
[http://www-nds.iaea.org/RIPL-2/]. 

[26] K. Govaert, F. Bauwens, J. Bryssinck, D. De Frenne, E. Jacobs, 
W. Mondelaers, L. Govor, and V. Y. Ponomarev, Phys. Rev. C 
bf57, 2229(1998). 

[27] A. J. Koning, S. Hilaire, and M. C. Duijvestijn, "TALYS-1.2", in 
Proceedings of the International Conference on Nuclear Data 
for Science and Technology, April 22-27, 2007, Nice, France. 
Editors: O. Bersillon, F. Gunsing, E. Bauge, R. Jacqmin, 
and S. Leray, EDP Sciences, 211 (2008). Available online at 
[http : //www. taly s . eu/] . 

[28] S. Goriely S. ffilaire, and A.J. Koning, Phys. Rev. C 78, 064307 
(2008). 

[29] A. J. Koning and J. P. Delaroche, Nucl. Phys. A 713, 231 

(2003). 

[30] P E. Koehler, R. R. Spencer, K. H. Guber, J. A. Harvey 
N. W. Hill, and R. R. Winters, in Proceedings of the Interna- 
tional Conference on Nuclear Data for Science and Technol- 
ogy, Trieste, Italy, 1997. Editors: G. Reffo, A. Ventura, and 
G. Grandi (ItaUan Physical Society, Bologna, Italy, 1997). 

[31] R. L. Macklin, T. Inada, and J. H. Gibbons, Nature (London) 
194, 1272 (1962). 

[32] J. Nishiyama, M. Igashira, T. Ohsaki, G. N. Kim, W. C. Chung, 
and T. I. Ro, J. Nucl. Sc. Tech. (Tokyo), 45, 352 (2008). 



13 



[33] K. Wisshak, F. Voss, Ch. Theis, F. Kappeler, K. Guber, L. Kaza- 

kov, N. Komilov, and G. Reifo, Phys. Rev. C 54, 1451 (1996). 
[34] V. M. Timokhov, M. V. Bokhovko, M. V. Isakov, L. E. Kaza- 

kov, V. N. Kononov, G. N. Manturov, E. D. Poletaev, and 

V. G. Pronyaev, Yad. Fiz. 50, 609 (1989). 
[35] E. J. Winhold, E. M. Bowey, D. B. Gayther, and B. H. Patrick, 

Physics Letters 32B, 7 (1970). 
[36] P Adrich et al, Phys. Rev. Lett. 95, 132501 (2005). 
[37] G. Tertychny, V. Tselyaev, S. Kamerdzhiev, F. Grilmmer, S. 

Krewald, J. Speth, A. Avdeenkov, and E. Litvinova, Phys. Lett. 

B 647 104 (2007). 
[38] N. Tsoneva and H. Lenske, Phys. Rev. C 77, 024321 (2008). 
[39] N. Paar, Rep. Prog. Phys. 70 (2007), p. 691-793. 
[40] I. Daoutidis (private communication). 

[41] J. Endres, E. Litvinova, D. Savran, P. A. Butler, M. N. Harakeh, 



S. Harissopulos, R.-D. Herzberg, R. Kriicken, A. Lagoyannis, 
N. Pietralla, V. Yu. Ponomarev, L. Popescu, P. Ring, M. Scheck, 
K. Sonnabend, V. I. Stoica, H. J. Wortche, and A. Zilges, Phys. 
Rev. Lett. 105, 212503 (2010). 

[42] P. Van Isacker, M. A. Nagarajan, and D. D. Warner, Phys. Rev. 
C45, R13 (1992). 

[43] D. Sarchi, P F. Bortignon, and G. Colo, Phys. Lett. B 601, 27 
(2004). 

[44] E. Litvinova, P. Ring, and V. Tselyaev, Phys. Rev. Lett. 105, 

022502 (2010). 
[45] A. Tonchev (private communication). 
[46] Klimkiewicz et al, Phys. Rev. C 76, 051603(R) (2007). 
[47] A. Tonchev et al, Phys. Rev. Lett. 104, 072501 (2010). 



